Origin of Gap Anisotropy in Spin Fluctuation Models of the Fe-pnictides 
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We discuss the large gap anisotropy found for the Aig (s-wave) state in RPA spin-fiuctuation 
and functional renormalization group calculations and show how the simple arguments leading to 
isotropic sign-switched s-wave states in these systems need to be supplemented by a consideration 
of pair scattering within Fermi surface sheets and between the individual electron sheets as well. In 
addition, accounting for the orbital makeup of the states on the Fermi surface is found to be crucial. 



PACS numbers: 

With one or two exceptions, the newly discovered Fe- 
pnictide superconductors are created by doping parent 
materials which manifest a magnetically ordered ground 
state. Since critical temperatures are as high as 56K, and 
experiments have shown enhanced magnetic response rel- 
ative to ab initio electronic structure calculations, it is 
natural to consider spin fluctuation type pairing mod- 
els for these systems. There have been several calcula- 
tions of this general type, which have reported gaps with 
significant anisotropy on the multisheeted Fermi surface 
of these materialsii^i^i^i^i^ii^i^i^. For the most part, a 
ground state with Alg symmetry has been identified, 
with a nearby d.j;2_y2 (in the unfolded Brillouin zone) 
pairing eigenvalue. These calculations involve fairly real- 
istic representations of the electronic structure near the 
Fermi surface, and are sufficiently complicated that the 
precise physical effects leading, e.g. to a particular pair- 
ing channel or to gap anisotropy can be obscured. 

At the same time, there is a simple argument origi- 
nally given by Mazin et al.^^ which suggests that the or- 
der parameter in the system should have a sign-switched 
s-wave structure^iiSiiiiii^iii^. It is based on the argument 
that the scattering of pairs between the a-hole Fermi 
surfaces around the F point and the /9-electron Fermi 
surfaces around the X(7r, 0) and y(0,7r) points of the 
unfolded Brillouin zone (see Fig. 1) is dominant in the 
system due to the near-nesting of the small sheets. In 
this picture, the system can maximize its condensation 
energy by forming isotropic order parameters but switch- 
ing sign between the a and (3 sheets to take advantage 
of the interband pair scattering. This change in sign 
has the added benefit of reducing the short range on site 
Coulomb repulsion. It is therefore at first sight sur- 
prising to find that RPA spin-fiuctuation calculationsiii^ 



as well as functional renormalization group studies^ find 
a highly anisotropic Aig s-wave gap. Here we investi- 
gate this, as well as explore the question of what the 
anisotropy can tell us about the pairing interaction. One 
answer that has been given is that the momentum de- 
pendence of the fiuctuation-exchange pairing interaction 
can drive the anisotropy^. However, we will see that the 
orbital make-up of the states on the Fermi surface and 
the suppression of the short range Coulomb interaction^ 
also play key roles in favoring an anisotropic gap. 

In the following we examine the pairing strengths for 
processes which involve a pair of electrons scattering on 
and between the four Fermi surfaces shown in Fig. [1] 
These were calculated using a 5-orbital {dxz, dyz, d-xy, 
dx2-y2, d3z2_r2) tight-binding fit to the DFT bandstruc- 
ture calculations of Cao et aU^. The tight-binding pa- 
rameters for an orbital basis that is aligned parallel to 
the nearest neighbor Fe-Fe direction are given in the ap- 
pendix of Ref. [l^. Here we will label the two hole Fermi 
surfaces ai and 0:2 and the two electron Fermi surfaces 
/3i and (32 as indicated in Fig. [1] 

The orbital weights a* (A:) of the states on the various 
Fermi surfaces are shown in Fig. [2] and illustrated by the 
colors in Fig.[TJ Here i designates the Fermi sheet ai, 0:2, 
and P2 and t the orbital (dxz, dyz,. . .). The dominant 
orbitals contributing to the ai and a2 sheets are the d^z 
and dyz orbitals. The upper and lower parts of the ai 
sheet are dyz-like and the left and right hand sides have 
dxz character. The opposite behavior is seen on the 0:2 
sheet. On the /3i sheet, the upper and lower surfaces have 
dominantly dyz character while the ends along the kx-s-'xis 
have dxy character. Similarly, the /32-sheet is made up of 
dxz along the sides and dxy on the ky-axis end o-^^'^^i^*" . 

Now consider the scattering of a pair (fc' |, —k' |) on 
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FIG. 1: (color online) The Fermi surface of the 5-orbital tight- 
binding modeli^. The main orbital contributions are shown 
by the following colors/symbols: d^z (red/solid circles), dyz 
(green/open circles), d^y (blue/diamonds). 
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the ai Fermi surface to a pair {k t, — fc i) on the f3i Fermi 
surface. The strength of this scattering depends upon 

Ti,{k.k')= (1) 
<f (-fc)<:*(fc)Re k\ 0)] al^ (fc')<, i-k') 

stpq 

and involves orbital weight factors and, within the fluc- 
tuation exchange approximation the orbital dependent 
vertex 
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Here, the momenta k and k' are restricted to the differ- 
ent Fermi surface sheets Ci with k G Ci and k' G Cj. 
The interaction matrices W and W contain the onsite 
intra- and inter-Coulomb interactions along with the ex- 
change couplings and x?^^ and x^^^ are the RPA spin 
and charge orbital susceptibilities. We will use numeri- 
cal results obtained from earlier work22 but our results 
are not dependent on the precise values of the param- 
eters. Rather, what is important to note is that the 
dominant pairing interaction is found to arise from the 
spin- fluctuation term ^U^x^^^U" and the short range 
Coulomb contact interactions jW^ and oppose the 
pairing. As noted by Mazin and Schmalian^'*, for a sim- 
ple two-band model with equal density of states, the short 
range Coulomb repulsion vanishes for a sign switched s- 
wave gap. However, as they point out, when there is an 
asymmetry in the density of states, some fraction of the 
Coulomb-repulsion remains. We will see that it is both 
the orbital weight factors and the further suppression of 



FIG. 2: (color online) The orbital weights as a function of 
the winding angle (j> on the different Fermi surface sheets^. 
The different colors/lines refer to d^z (red/solid lines), dyz 
(green/dashed lines), d^y (blue/dash-dotted line), d^2_y2 
(yellow/short-dashed line). 



the Coulomb interaction that lead to the anisotropy of 
the gap in the realistic case. 

First we discuss the effect of the orbital matrix ele- 
ments. The scattering strength ry(fco,fc) for two dif- 
ferent pairs on the ai-sheet with ko = (^^,0) and 
ko = (0, kp) as a function of momentum k are illustrated 
in Fig. [3l For the former. Fig. [3^, the dominant scat- 
tering is to dxz{k I, —k J.)-pairs on the (32 Fermi surface 
and for the latter. Fig. to dy^-pairs on the /3i Fermi 
surface. Note that the scattering strength is not simply 
a consequence of nesting; instead, it reflects the orbital 
weight structure factors. The d^zik' f, —k' |) pairs on ai 
scatter more strongly to d^zik t, —k i) pair states on (32, 
than to pair states on the Pi Fermi surface which involve 
other orbitals. This means that while it is in general fa- 
vorable to have a sign change between the gaps on the 
ai and P2 Fermi surfaces, the essential thing is to have 
a sign change of the gap between the red (dxz) regions 
of the «! Fermi surface and the red {dxz) parts of the 
(32 Fermi surface shown in Fig. [TJ Likewise, one needs 
a sign change between the green (dyz) regions of the ai 
Fermi surfaces and the green (dyz) Pi Fermi surface. It is 
not important to maintain this sign change in the yellow 
{dxy) regions of the P Fermi surfaces. We will in fact 
see that the magnitude of the gap is larger on the dyz 
portion of the Pi Fermi surface and smaller on the dxy 
parts. This actually leads to a small increase in the a~Pi 
pairing compared to the isotropic sign-switched gap. 
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FIG. 3: (color online) The strength Tij{ko, k) associated with 
scattering a pair from the ai Fermi surface with momenta a) 
ko = (kpi, —kpi) and b) ko — {kpj, —kFj) as indicated by 
the black circles as a function of momentum k. The orbital 
weight factors favor scattering from dxz to dxz and dy^ to dy^ 
orbital states. 



In addition, as we will discuss, if there are other inter- 
actions such as the local Coulomb interaction or scatter- 
ing processes such as inter-Fermi surface P1-P2 scatter- 
ing, reducing the magnitude or even changing the sign of 
the gap on the dxy part of the (3 Fermi surfaces can lead 
to an additional reduction of the short range Coulomb 
interaction and an enhancement of the pairing. In order 
to explore this latter effect, we need to obtain a more 
detailed accounting of the various contributions to the 
pairing strength. As discussed in Ref. [isl . for a given gap 
function g{k), the effective pairing strength is determined 



A[.9(fc)] 
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Here vp{k) = \VkEi{k)\ for fc on a given Fermi surface 
Ci. In the following, we will normalize the gap function 
g{k) such that the denominator is equal to the total one- 
electron density of states, 
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Then we can decompose A into its contributions from the 
different inter- and intra-Fermi surface scattering pro- 
cesses. 
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First consider for two limiting cases (1) the opti- 
mal RPA g{k) found as the variational solution of Eq. [3] 
and (2) a sign-switched s-wave with g^-^ = ga^ — 1 and 
.9/3i = 3& = —1- The results are shown in Figs. [Hand 
The total pairing strength is significantly larger for the 
RPA g{k). From the breakdown of the various contri- 
butions one sees that while the total a-/? contribution is 
larger, 0.23 compared with 0.21, for the sign-switched gap 
function, the larger negative contribution of the Fermi 
surface A's as well as the larger negative /3i-/32 contri- 
bution overcome this and A for the sign-switched gap is 
considerably smaller than the optimal A. In Fig. [H one 
can see that the anisotropy of the RPA solution is such 
that there are nodes on the /^-sheets^- . As we will discuss, 
whether this happens or not depends upon the parame- 
ters. The interplay of the orbital weights in the pairing 
interaction, Eq. [3l and the reduction of the intra (3i-(3i 
Coulomb interaction, as well as the inter /3i-/?2 scattering 
lead to the anisotropy. 

To illustrate this, consider the simple parameterization 



of an anisotropic gap with g^^ 



1 and 



5/3 



-a(l - r cos 26*) 



(7) 



1/2 

Here a = (2/(2 + r^)) so that the normalization con- 
dition, Eq. (dl is satisfied. When r = 0, gp — ~1 and we 
have the sign-reversed state previously discussed. Then 
as r increases, the gap function becomes anisotropic on 
the P Fermi surfaces as shown in Fig. [HI In Fig. [7] we plot 
the total pairing strength A and some of its components 
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FIG. 4: (color online) The RPA gap function a) plotted on 
the Fermi surfaces (red/solid circles positive and blue/open 
circles negative), (b) g(k) as a function of angle and (c) bar 
graph of the various interaction components Xij. Note that 
the off-diagonal terms contribute to the total A, Eq. ([5]) with 
weights of two or four as indicated. 



FIG. 5: (color online) The sign-switched gap function (similar 
to Fig.g. 



versus r. Initially, as r increases, the gap on the f3 Fermi 
surface sheets becomes anisotropic, and the total pairing 
strength A increases. The slight increase in Xa-p reflects 
the increase in the amplitude of the gap in the regions 
where the d^z and dy^ orbital weights are largest. The 
additional increase in A arises from the suppression of 
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FIG. 6: (color online) The phenomenological anisotropic gap 
function on the /3 Fermi surface versus angle for several dif- 
ferent values of r. 



the intra- and inter-Coulomb repulsion associated with 
the /?i and /52 Fermi surfaces due to the gap anisotropy. 
Finally, as r increases further, the reduction of the a-fi 
pairing contribution due to the anisotropy becomes larger 
than the suppression of the Coulomb interaction and the 
total pairing strength A decreases. For the interaction 
parameter set that we are using, this occurs for r ~ 1.5, 
so that there are well developed nodes on the /3-Fermi 
surfaces. However, if the intra- and inter-scattering on 
the /3-Fermi surfaces were reduced, the nodes could be 
lifted. 



Conclusion 

The anisotropy of the A\g gap on the electron Fermi 
surfaces found in RPA and numerical functional renor- 
malization group studies has been shown to arise from 
an interplay of three sources: (1) the variation of the 
weighting of the different d-orbitals on the Fermi surfaces, 
(2) the need to suppress the Coulomb repulsion and (3) 
the need to reduce the effects of the repulsive scatter- 
ing between the electron-/? sheets. The orbital weight 
variation is familiar in other multi-orbital superconduc- 
tors such as MgBai^. In essence, the pairing interac- 
tion is strongest between fermions in near neighbor dyz 
orbitals along the x-direction and near-neighbor d^^ or- 
bitals along the y-direction. The pairing associated with 
the dxy near neighbor orbitals is weaker. At the same 
time, the anisotropy leads to a reduction of the repulsive 



FIG. 7: (color online) The pairing strength and its com- 
ponents versus r for the anisotropic gap shown in Fig. |6] 
Here A is the total pairing strength. The intra- and inter- 
OL\ and 02 Fermi surface contributions are given by = 
Aqiqi + Aa2ci2 + 2Aaici2 thc lutra-/? Fermi surface contri- 
bution by A^ = A^j^j -I- A^2i32 ■ The contribution to the pairing 
comes from Aq-^ = '^XaxPx + 2Ac,2/3i , while for this A\g gap 
the inter-/3i — /32 contribution opposes the pairing for small r 
and is neutralized for a large anisotropy. The red dashed line 
indicates the optimal degree of anisotropy, r = 1.5. 



Coulomb interactions and the inter-Fermi surface /3i-/32 
scattering. As discussed, there is a balance between these 
effects which determine whether the anisotropy is suffi- 
ciently large that there are nodes on the /3-sheets. The 
fact that the dominant orbital weights on the a\ and 
Q!2 Fermi surfaces are associated with the d^z and dyz 
orbitals leads to a more isotropic gap on these sheets. 

The results that we have discussed are based upon a 
weak coupling spin-fluctuation approach. In particular, 
the various calculations have used an RPA form for the 
pairing vertex. Here the bandstructurc and filling, along 
with the relative strengths of the onsite Coulomb and ex- 
change interactions enter. A key feature of this approach 
is that the gap exhibits large anisotropics and that nodes 
may appear on the /3-Fermi surfaces in the s-wave A\g 
state. However, one can imagine that variations in the 
parameters can alter the degree of anisotropy. An impor- 
tant implication of this is that it may provide an explana- 
tion for the wide variety of experimental indications that 
nodes are present in some cases and not in others. An- 
other explanation for these apparent discrepancies could 
be the degree of disorder in different samples as discussed 
in Refs. ftolilll^iil. 
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